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ABSTRACT 


Liquid  metal  combustion  has  been  used  as  the  primary  energy  source  for 
underwater  vehicle  propulsion  because  of  its  energetic  and  closed-cycle 
characteristics.  This  can  be  achieved  by  injecting  gaseous  oxidizer  (SFg)  in  a 
molten  metal  (Li)  bath.  The  present  paper  describes  a  new  multi-fluid  (MF) 
combustion  model  for  multiphase  submerged  diffusion  flames.  Each  fluid  has 
separate  transport  equations  which  are  linked  to  other  fluids'  transport 
equations  by  additional  source  terms.  The  present  MF  model  results  for  a 
submerged  Li(f)-SFg(g)  flame  gave  a  much  longer  plume  length  compared  to  the 
corresponding  local ly-homogeneous  flow  (LHF)  model  results  and  is  in  good 
agreement  with  the  experimentally  measured  plume  length. 
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1.  INTRODUCTION 


Stored  chemical  energy  propulsion  systems  (SCEPS)  using  Li/SFg  reactants 
have  recently  received  a  great  deal  of  attention  as  a  means  of  producing  high 
energy  per  unit  mass  of  reactant  [1-9].  The  "SCEPS"  system  utilizes  a 
liquid-metal  fuel  combustor  as  part  of  a  closed  steam-power  cycle.  The  combustor 
uses  a  liquid  metal  bath  (fuel)  through  which  a  high-momentum  gas  jet  (oxidizer) 
is  injected  forming  a  reacting  multiphase  submerged  jet.  Lithium  (Li)  has  been 
used  as  the  liquid  metal  while  sulfur  hexafluoride  (SFg)  gas  or  other  halogen  gas 
is  used  as  the  oxidant.  In  this  case,  flame  temperatures  as  high  as  4000  K  may 
be  encountered  making  detailed  flame  structure  measurements  currently  very 
difficult,  if  not  practically  impossible. 

Computer  simulation  of  liquid  metal  fuel  combustors  has  recently  been 
adopted  as  the  next  possible  technique  for  obtaining  detailed  information 
regarding  the  multiphase  submerged  flame  structure  [1-11].  All  the  current 
computational  models  use  the  LHF  model  as  a  basis  for  predicting  the  submerged 
diffusion  flame  with  the  obvious  assumption  of  infinitely  fast  transport 
mechanisms  [phase  equilibrium  and  same  velocity  and  temperature  of  all  phases, 
etc.]  between  the  reacting  phases.  This  very  assumption  leads  to  a  substantially 
underpredicted  plume  length  for  the  Li-SFg  submerged  flame  [9]. 

Presently,  other  than  the  LHF  model,  very  few  combustion  models  can  be  used 
to  predict  the  turbulence  structure  of  multiphase  submerged  reacting  plumes. 
Practically  all  the  contemporary  non-submerged  models  (stochastic  separated- flow 
model,  etc.)  developed  for  liquid-fuel  spray  flames  are  unfortunately 
inapplicable  to  the  submerged  reacting  plumes  because  of  the  flow  regime  change 
and  condensation  of  products  and  excess  evaporated  fuel.  When  a  gaseous  oxidant 
is  injected  into  a  liquid  fuel,  the  flow  regime  begins  with  a  single-phase, 
gaseous  flow  pattern  at  the  nozzle  exit.  As  the  liquid  fuel  is  entrained  to 
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react  with  the  oxidant,  a  liquid  droplet  flow  pattern  may  evolve  first.  However, 
the  mixture  temperature  is  still  low  enough  that  the  product  may  appear  in  a 
solid  phase  [9]  such  that  the  droplet  flow  regime  may  be  composed  of  gas,  liquid 
and  solid  phases.  Further  downstream,  with  sufficient  entrainment  and  reaction, 
a  gaseous  flow  pattern  may  reappear,  which  may  be  followed  by  another  droplet 
flow  and  then  a  bubbly  flow  pattern  due  to  condensation  of  vapor  products  and 
excessively  evaporated  fuel  vapor.  Eventually  it  should  turn  into  condensed 
phases  at  the  tip  of  the  plume.  Though  the  above  description  of  flow  patterns 
may  be  rather  speculative,  the  submerged  reacting  jets  are  therefore  distinctly 
different  from  the  better  known  two-phase  spray  jets,  in  which  droplets  from  the 
spray  nozzle  simply  evaporate  and  react  to  form  non-condensable  gaseous  products 
with  no  subsequent  condensation.  It  is  somewhat  similar  to  the  conventional  two- 
phase  flow,  see  for  example  Ishii  [12],  in  which  the  void  fraction  can  range 
between  zero  and  one,  flow  regimes  can  change  and  the  interaction  between  two 
phases  is  strong.  Consequently,  the  continuum  concept  of  the  two-fluid  model  in 
the  field  of  two-phase  flow  is  adopted.  However,  it  is  necessarily  extended  to 
reacting  multiphase  flow  to  account  for  not  only  velocity  slip  between  phases  but 
also  phase  nonequilibrium  between  phases  and  turbulence-reaction  interaction 
through  a  PDF  (probability  density  function)  approach.  The  results  are  compared 
with  that  of  the  prevailing  LHF  model  and  available  experimental  data  to  validate 
the  proposed  model  for  multiphase  submerged  combustion.. 

2.  THE  PHYSICAL  MODEL 

The  proposed  physical  model  treats  the  multiphase  submerged  flame  as  a 
number  of  "fluids"  or  phases  (i.e.,  fluid  1,  fluid  2,  ...,  fluid  k,  ...,  fluid 
n).  The  total  number  of  fluids  is  n;  any  flow  or  thermochemical  property  (0) 
pertaining  to  fluid  k  (1  <  k  <  n)  is  designated  by  a  "k"  superscript,  e.g.,  0*^). 
The  last  two  fluids  (fluid  (n-1)  and  fluid  n)  are  respectively  dispersed 
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condensed  combustion  products  (such  as  LiF(l),  Li2S(f)  and  LiF(s))  and  a 
continuous  gas  phase  (such  as  SFg(g),  Li(g),  Li2F2(g),  Li3F3(g)  and  Li2S(g), 
etc.)*  The  remaining  fluids  (1  <  k  <  (n-2))  are  dispersed  Li  droplet  phases; 
each  fluid  is  characterized  by  a  mean  droplet  diameter  with  the  smallest 
diameter  at  k  >  1.  For  compactness,  all  the  modeled  time-averaged,  steady-state, 
conservation  equations  given  below  are  expressed  in  compact  tensor  notation 
although  they  are  actually  solved  in  their  cylindrical  polar  form.  In  the 
governing  equations,  terms  involving  p'  were  ignored.  In  the  I’mit  of  phase 
equilibrium  and  non-slip  between  phases,  the  equations  reduce  to  the  Favre- 
averaged  form  of  the  LHF  approximation  [8]. 

a)  The  mass  conservation  equation  for  fluid  k: 


-  C  (1) 

where  the  comma  suffix  denotes  differentiation  with  respect  to  the  spatial 
coordinates,  is  the  turbulent  viscosity  of  fluid  k;  0*^^  is  a  turbulent  Schmidt 
number;  5*^  is  the  mean  volume  fraction  of  fluid  k;  u^.  is  the  mean  velocity  in 
the  i -direction;  is  the  density;  is  a  volumetric  mass  source  term  because 
of  mass  transfer  between  the  fluids  due  to  evaporation,  condensation  and  droplet 
inter-fluid  transfer  as  droplets  decrease  in  size.  is  given  as. 


-  m  -  (m  (1  <  k  <  (n  -  2)) 

and  for  k  -  n  (carrier  fluid) 

'  n-1 

-  2  m^^a*  -  m  (k  »  n) 


(2) 


(3) 
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In  the  above  is  the  volumetric  liquid  phase  mass  transfer  from  fluid  (k+1) 
to  fluid  k.  It  can  represent,  for  example,  the  rate  of  mass  of  droplets  in  size 
(k+1)  group  that  is  shifted  to  size  (k)  group  due  to  evaporation.  Similarly,  m 
is  the  mass  flowing  out  (as  liquid  phase)  from  fluid  k  and  m^  is  the 
evaporation  rate.  For  the  condensing  combustion  products  (k=n-l)  Equation  (1) 
is  actually  used  to  compute  the  condensation  rate  source  term  as  will  be 

explained  later  in  Section  2.5.  The  last  term  in  Equation  (3)  represents  the 
condensation  rate  which  is  a  loss  of  mass  from  the  carrier  phase  to  the  (n-1) 
phase,  assuming  that  the  reaction  has  taken  place  in  the  gas  phase  only. 

b)  Modeled  momentum  equation: 

^  -k„  /„kw-k  .  -k  2  -k  ,  . 

y  >  u,Uj  -  a  (Hf  to  v^/0  )(u,_j  t  Uj_,  -  3 

1  ;k^kjt  :;k,k/-k-k  .  -k-k  ,  .  ki  -k- 

-  3  0  k  «ij  -  0  “  "J 

*  fSj  *  (4) 


where,  k*^  is  the  turbulence  kinetic  energy;  and  are  the 
kinematic  turbulent  viscosities  of  the  carrier  phase  (k=n)  and  fluid  k;  the  mean 
drag  and  mass  transfer  forces  and  may  be  given  as  follows  for  1  <  k  < 
n-1  (dispersed  fluids) 
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and  for  k  »  n  (carrier  fluid) 


n-1 

2 
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(8) 


Finally,  the  global  continuity  must  be  included 

2  =  1  (9) 

k 

In  the  above  equations,  (a  is  a  turbulent  Schmidt  number;  is  the 

molecular  viscosity;  is  the  gravitational  acceleration  in  the  j-direction; 
is  a  reference  density;  p  is  the  pressure;  f'^p  is  a  drag  function  which  may  be 
written  as  [13] 

fp  -  18  Z  m5/(d'^)^  (1  <  k  <  n)  (10) 

where  is  a  mean  droplet  size  of  "fluid  k";  Mf"  is  the  carrier  phase  molecular 
viscosity  (i.e.,  for  k  =  n);  Z  is  a  function  of  the  droplet  Reynolds  number  Re'^ 
[13]  and  the  droplet  transfer  number  B  [14] 

Z  =  [1  +  0.15  (Re*^)®-®®^  +  0.02  ReV(l  +  42500/(Re'^)^-^®]/ 

(1  +  (11) 

where,  Re*^  a  p*’D*^|u!)  -  u*[  |//i,;  is  the  carrier  phase  density. 

The  turbulent  Schmidt  number,  <^y  reflects  the  coupling  effect  between  the 
gas  and  liquid  momentum  fields.  It  can  be  calculated  as  [15] 

.  1  +  I  jVtc^kfJd+iV/A")]  (k^ti)  (12) 

u 

and  for  k  -  n,  a  *  1  by  definition. 
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In  the  transport  equations,  only  second  order  correlation  terms  (e.g., 


a a ‘‘u/")  were  retained  [16]. 

2.1  The  Turbulence  Model 

A  multi -fluid  k-6  model  is  adopted  where  transport  equations  for  the 

kinetic  energy  of  turbulence  k  and  its  dissipation  rate  e,  of  the  carrier  fluid, 

are  solved  [17],  namely 

k  i]^i  =  a"  (G  +  AG) 

-  A"  (c  +  Ac)  (13) 

A"u-c  .j  -  [(/ita%f)  Cj] j  =  a"  (f)  Cj  (G  +  AG) 

-  A"  (f)(C2€  +  C3AC)  (14) 

where,  G,  AG  and  Ac  are  the  single  phase  generation  term,  extra  generation  and 
dissipation  terms  due  to  volume  fraction  and  velocity  fluctuations  [17].  of 
the  carrier  phase  is  defined  as 

l-t  »  (15) 

where  C^,  Cj,  C^,  C3,  and  are  constants  of  the  model  [17]. 

2.2  Inter-Fluid  Mass  Exchange  Model 

The  different  fluids  flowing  through  the  submerged  plume  exchange  mass, 
momentum  and  energy.  The  exchange  of  mass  is  mainly  controlled  by 
evaporation/condensation  inside  the  plume. 

Droplets  are  formed  at  the  interface  between  the  injected  gaseous  oxidant 
(SFg)  and  the  liquid  metal  bath  (Li).  The  size  of  the  formed  droplets  is 
dependent  on  the  critical  wavelength  of  the  disturbed  gas-liquid  interface  [18]. 
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This  wavelength  is  a  function  of  the  nozzle  jet  velocity  and  the  liquid  phase 
viscosity  and  surface  tension  [18].  The  entrained  droplets  then  decrease  in  size 
as  a  result  of  evaporation.  Hence,  they  flow  (in  droplet-size  space)  from  a 
fluid  of  a  larger  droplet  size  to  one  with  smaller  droplet  size  (i.e.,  from  fluid 
(k+1)  to  fluid  k) . 

The  conservation  equations  (1  and  4)  for  fluid  k  are  affected  by  the 
evaporation  rate  m^ outflow  of  mass  m  and  the  inflow  of  mass  m  due  to  the 
decrease  in  droplet  diameter  (doVdt).  The  expressions  for  lii^  and  m for  k 
<  n-2,  are  [19] 


"7  ■  ^  <e^V[2(D'‘)^]  (16) 

.  (d'‘  -  40^2)^  Kj^'‘/[2(d'')^5d'']  (17) 


where,  S&  is  the  size  increment  for  fluid  k  and  is  the  forced  convection 
evaporation  constant.  It  can  be  calculated  as  [14] 


.  [1  ^  at  (1  ,  B) 

®  (1  +  r 


(18) 


where,  Pr  is  the  gas  film  Prandtl  number;  T  (s  Hf/Pr)  is  an  exchange  coefficient; 
the  droplet  transfer  number  B  is  Cp  (T  -  Tjj'^)/l'^;  Cp  is  the  gas  film 
constant-pressure  specific  heat;  and  T  are  the  boiling  temperatures  of  fluid 
k  and  the  gas  (fluid  n)  temperature;  is  the  latent  heat  of  vaporization. 
Equation  (17)  can  also  be  used  for  the  inflow  of  droplets  from  k  +  1  to  k 

fluid,  if  k  is  replaced  by  k  +  1. 

As  the  temperature  of  combustion  products  decreases  at  the  end  of  the 
plume,  new  droplets  can  be  formed  as  a  result  of  condensation.  Therefore,  mass 
is  transferred  from  fluid  n  to  fluid  n  -  1.  It  is  not  expected  that  the 
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condensed  product  will  evaporate  again  hence  lii^  =  th  =  0.  When  k  =  n,  rh^ " 
=  0  (gas  phase)  while  lii "  is  interpreted  as  the  condensation  rate  per  unit  volume 
of  the  gas  phase.  The  latter  can  be  calculated  as  explained  in  section  2.5. 

Phase  Non-Eoui librium  Variable 

Due  to  the  finite  evaporation/condensation  rate,  the  mean  mass  fraction  of 
the  liquid  fuel  (Li),  ^  (s  ^5*'),  will  not  be  in  equilibrium  with  the 

k 

corresponding  gas  phase.  Therefore,  a  mean  non-equilibrium  variable  e  may  be 
defined  as 

n-2 

i  -  t(Y„  -  V^)/(Y„  -  Y^jJ^i.  (19) 

where  {Yy)LL  and  are  the  nonreacting  and  reacting  (equilibrium)  liquid  Li 

mass  fractions.  is  to  be  obtained  from  a  full  equilibrium  Li-SFg  state 

relationship  (e  »  1)  while  (yjii  is  simply  (1  -  f),  f  being  the  mean  mixture 
fraction. 


2.3  The  Combustion  Hodel 

In  the  present  combustion  model,  chemistry  is  assumed  fast  enough  so  that 
the  multiphase  reacting  mixture  is  in  chemical  equilibrium.  However,  the 
evaporation/condensation  rate  is  not  so  fast  for  phase  equilibrium  to  prevail  at 
all  points  inside  the  plume.  Therefore,  the  instantaneous  state  of  the  reacting 
mixture  is  a  function  of  f  and  the  non-equilibrium  variable  e  (Eq.  19). 

The  time-averaged  value  of  any  thermophysical  property  (^)  solely  dependent 
on  f  and  e  may  be  computed  as  [19] 


}  }  ^(f.e)  P(f)  P(e)  df  de 
0  0 


(20) 
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where  P(e)  is  determined  from  e  as 

P(e)  =  (1  -  e)  5(e)  +  e  5(1  -  e)  (21) 

and  P(f)  is  the  probability  density  function  characterized  by  f  and  g  (s  f'f'). 
The  above  two  equations  yield 

1 

^  =  J  [(1  -  i)  ^(f,o)  +  i  0(f,l)]  P(f)  df  (22) 

0 

where  0(f,o)  is  the  value  of  0  with  no  evaporation  (e  =  o)  while  0(f,l)  is 
related  to  f  through  full  thermodynamic  equilibrium  state  relationships  (e  =  1). 
A  beta- function  pdf  is  adopted,  namely 

P(f)  =  f®'^  (1  -  f)^'V  J  (1  -  f)^'^  df  (23) 

0 


where  a  and  $  are  given  as 
a  =  f  [f  (1  -  f)/g  -  1] 

=  a  (1  -  f)/f 

f  and  g  are  to  be  computed  from  their  transport  equations. 


(24) 

(25) 
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An  exact  instantaneous  two-fluid  mixture  fraction  transport  equation  was 
introduced  by  Chan  and  Abou-Ellail  [16].  This  equation  can  be  extended  to  the 
multi -fluid  case.  The  modeled  MF  form  reads 


[p  u.f 


^^f,eff^  ^,i^i 


where 


2  n)/S  +  (?  +  ”2^ 
k*l  ^  ^  k-1  ^  ^ 


(26) 
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where  p  (a  2  is  the  mixture  density;  u.  (a  2  is  the  mixture 

k  k  _ 

velocity;  is  an  effective  exchange  coefficient;  is  an  additional 

source/sink  term  [19]  resulting  from  the  interaction  of  the  multi-fluids. 
Similarly,  the  g  transport  equation  can  be  written  as 


“l9  ■  Cg.eff)  9.l5.1 


Si  "t 


-  Cg2M</k 


where  r 

g> 


eff  “ 


2  n)/S  +  (^  +  a^/a'^)nJa 
k=l  ^  ^  k=l  ^  ^ 


(27) 


In  the  above  equations,  o^,  a^,  C^j  and  Cg2  constants  of  the  physical  model 
while  AG^  is  an  extra  generation/dissipation  term  [16,19]. 

The  model  constants  used  in  the  present  paper  are:  =  0.09,  Cj  =  1.44, 

C2  =  1.84,  C3  =  1.2,  Cgj  *  2.8,  Cg2  •  1.84,  “  0:7 ^  0^  *  1.0,  S^.  =  0.7, 

*  1.3  and  *  1.0. 

2.5  State  Relationships 

The  present  MF  reacting  model  utilizes  an  instantaneous  state  relationship 
which  is  a  function  of  the  mixture  fraction  f  and  the  phase  non-equilibrium 
variable  e.  A  full  equilibrium  (e  -  1.0)  state  relationship  for  Li(f)-SF5(g)  can 
be  constructed.  A  complete  description  of  the  equilibrium  computational 
algorithm,  which  is  mainly  based  on  the  Gibbs  free-energy  minimization  approach, 
is  given  by  Chan  et  al.  [8,9,20]. 

For  the  no  evaporation  case  (e  -  0),  the  liquid  and  gas  mass  fractions  and 
the  mixture  enthalpy  are  linearly  related  to  f  (i.e.  0(f,O)  -  ^(0,0) (1-f)  + 
^(l,0)f),  as  in  this  case  the  instantaneous  state  simply  represents  a  mixing 
process  with  no  chemical  reaction.  From  the  mixture  enthalpy,  mass  fractions  of 
the  liquid  and  gas  phases,  the  mixture  temperature  and  the  gas  and  liquid  volume 
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fractions  are  then  calculated  at  each  value  of  f  for  the  e  *  0  plane.  Constant 
e  planes  (o  <  e  <  1.0)  can  thus  be  constructed  from  the  e  =  1  and  e  »  0  planes 
(e.g.,  0(f,e)  =  (l-e)  ^{f,0)  +  e  0{f,l)). 

Figure  1  depicts  the  overall  equilibrium  Li{i)-SF5{g)  state  relationships, 
0(f,l),  at  atmospheric  pressure  and  a  liquid  lithium  temperature  of  1273  K  while 
the  injected  SFg  gas  temperature  is  298  K.  It  should  be  mentioned  here  that  the 
mixture  density  and  void  fraction  (gas-phase  volume  fraction)  values  given  in 
Fig.  1  are  only  used  in  the  LHF  model  calculations.  A  detailed  account  of  all 
species  concentrations  is  given  by  Chan  et  al.  [9].  Since  the  mean 
non-equilibrium  state  is  fixed  by  f  and  e,  of  the  condensed  products  can  be 
calculated  as  follows.  First  Eqs.  (19)  and  (22)  can  be  used  to  find  5"  of  the 
carrier  phase 

5"  .  }  [(1-i)  a"  (f,0)  +  e  a"  (f,l)]  P(f)df  (28) 

0 

where  af"(f,0)  and  o"(f,l)  correspond  to  the  e-0  and  e-1  planes  of  the  states 
relationships.  Since  all  volume  fractions  of  the  liquid  Li  droplets,  ^^(k-l  to 
n-2),  are  to  be  calculated  from  Eq.  (1),  then  Eq.  (9)  can  be  used  to  calculate 

-  1  -  V  -  ? 

k-1 

With  the  value  of  the  condensation  rate  source  term  S  and  hence 

m""(»s’^"'^Va")  can  be  computed  from  Eq.  (1)  for  k  -  n-1. 

3.  BOUNDARY  CONDITIONS 

The  present  predictions  are  limited  to  the  case  of  a  vertical  gaseous  SFg 
submerged  Jet  (at  298  K)  injected  into  molten  Li  (at  1273  K).  The  nozzle 
diameter  (d)  is  2  mm  while  the  inlet  gas  velocity  (u^)  is  325  m/s.  The  initial 

f 
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turbulence  intensity  is  taken  as  5%  of  the  mean  gas  velocity  while  the  turbulent 
length  scale  is  taken  as  3%  of  the  injector  radius.  The  surrounding  molten  Li 
is  assumed  stagnant.  The  entrained  mean  droplet  diameter  is  estimated  [18]  as 
60  /xm.  The  above  data  is  essentially  equivalent  to  the  experimental  Li(f)-SFg(g) 
combustor  data  of  Parnell  et  al.  [21].  The  number  of  fluids  (n)  is  taken  as  6; 
the  first  four  fluids  are  dispersed  pure  liquid  Li  with  mean  droplet  diameters 
of  10,  30,  50,  60  /un.  Fluid  5  is  occupied  by  the  condensing  combustion  products 
(70  im  mean  droplet  diameter)  while  fluid  6  is  only  for  the  gaseous  phase.  At 
the  nozzle  exit  plane  (k  <  6)  =  0  while  5*^  (k  =  6)  =  1.0.  At  the  free  stream 
boundary  (gas-liquid  interface)  »  O!^  »  ot^  »  a®  »  0.0  while  a*  »  1.0,  since 
fluid  4  is  the  only  fluid  available  for  the  60  /un  entrained  Li  droplets. 

In  addition  to  the  multi-fluid  (MF)  results,  similar  computations  were 
performed  for  the  LHF  model  to  be  used  as  a  basis  for  comparison.  In  this  case, 
at  the  nozzle  exit  plane,  all  volume  fractions  were  set  equal  to  zero  except  for 
the  carrier  phase  where  the  volume  fraction  was  set  equal  to  1.0;  the  density  is 
taken  as  that  of  the  mixture. 

4.  SOLUTION  PROCEDURE 

The  solution  procedure  is  based  mainly  on  an  extension  of  a  "multi -fluid 
SIMPLE”  algorithm  utilizing  an  iterative-marching  integration  technique 
[13,16,22].  The  main  variables  are  u*^,,  p',  k,  e,  f,  g  and  where  p'  is  a 
"multi-fluid  pressure  correction."  The  submerged  jet  is  overlaid  with  a 
non-uniform  2D  axisymmetric  grid  (700  axial  x  35  radial  nodes).  The  distances 
between  the  radial  nodes  are  varied  in  proportion  to  the  submerged  jet  spreading 
rate. 
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5.  RESULTS  AND  DISCUSSION 

The  predicted  axial  mean  profiles  of  the  centerline  mixture  fraction  and 
normalized  axial  velocity  for  the  MF  and  LHF  models  are  depicted  in  Fig.  2.  In 
Fig.  2,  u,  Ug,  0*^  and  are  axial  velocities  of  the  mixture,  gas  phase,  fluid 
4  and  SF„  at  nozzle  exit.  The  MF  model  results  show  that  both  f  and  u„  started 
to  decrease  much  further  downstream  (x/d  «  7)  than  the  corresponding  f  and  Ug  for 
the  LHF  model  (x/d  =1.5).  This  means  that  the  MF  model  predicts  a  much  longer 
plume  which  is  consistent  with  the  experimental  data  of  Parnell  et  al .  [21].  It 
is  interesting  to  notice  that  the  fluid  with  60  im  mean  droplet  diameter  (pure 
liquid  Li)  has  a  peak  axial  velocity,  near  x/d  »  6,  which  is  less  than  the 
prevailing  gas  velocity.  However,  for  x/d  >  70,  u*'  is  larger  than  Ug. 

The  existence  of  u*‘  in  the  whole  range  of  x/d  shown  in  Fig.  2  implies  that 
pure  liquid  Li  can  penetrate  to  the  centerline  of  the  plume  due  to  finite 
evaporation  rate.  Thus,  droplets  can  exist  at  the  centerline  even  though  their 
volume  fraction  e?  can  be  very  small  in  the  hot  region  of  the  plume. 

Figure  3  shows  the  axial  profiles  of  the  centerline  mean  gas  temperature 
T  and  the  mean  volume  (void)  fraction  of  the  gas  phase  for  the  MF  and  LHF  models. 
While  the  MF  model  predicted  a  peak  temperature  at  x/d  =  24,  the  LHF  model  peak 
tenqserature  occurs  at  x/d  =  1.5.  The  total  plume  length,  which  includes 
reaction,  evaporation  and  condensation  zones,  is  taken  at  ig  -  0.1.  The  present 
MF  model  predicts  a  plume  length  (L)  of  approximately  55  nozzle  diameters  while 
(L/d)  is  substantially  underpredicted,  by  the  LHF  model,  to  a  value  of  6.  The 
Parnell  et  al.  data  [21]  indicated  that  (L/d)  for  the  Li(f)-SFg(g)  system  is 
approximately  50  ±  10  depending  on  the  bath  temperature,  %  utilization  and  the 
submerged  Jet  orientation  with  respect  to  the  horizontal  direction.  Since  the 
MF  model  predicts  a  longer  plume  compared  to  that  of  the  LHF  model,  the  gradient 
of  (f,  J  Is  expected  to  be  smaller  and  therefore,  g  calculated  from  Eq.  (27)  will 
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also  be  smaller  resulting  in  a  narrower  pdf  profile.  This  explains  the  higher 
maximum  temperature  predicted  by  the  MF  model  as  can  be  seen  in  Fig.  3. 

Finally,  the  predicted  radial  profiles,  using  the  MF  model,  of  the  mean 
temperature  and  mean  mixture  fraction  are  depicted  in  Fig.  4.  for  x/d=20.  It  can 
be  seen  from  Fig.  4  that  the  maximum  temperature  (~3725K)  has  not  reached  the 
center  of  the  plume  yet  where  the  temperature  is  still  at  2100  K,  i.e.,  the 
reaction  zone  is  still  longer.  At  x/d»20  the  plume  outer  radius  occurs  at 
r/x~0.18,  based  on  f/f^  radial  profile  shown  in  Fig.  4. 

6.  CONCLUDING  REMARKS 

The  present  multi -fluid  (MF)  reacting  model  was  used  to  predict  the 
thermophysical  and  hydrodynamics  of  an  SFg  jet  submerged  in  molten  Li  fuel.  The 
predictions  for  the  plume  length  is  in  better  agreement  with  the  corresponding 
experimental  data,  which  could  be  considered  as  a  substantial  improvement  over 
the  local ly-homogeneous  flow  (LHF)  model.  The  latter  is  the  only  model  currently 
available  for  the  investigation  of  multiphase  submerged  combustion. 


-  14  - 


REFERENCES 


1.  HUGHES,  T.G.,  SMITH,  R.G.  AND  KILEY,  D.H.,  "Stored  Chemical  Energy 
Propulsion  System  for  Underwater  Applications,"  J.  of  Energy.  7,  pp.  128- 
133,  1983. 

2.  COOK,  L.P.  AND  PLANTE,  E.R.,  "Survey  of  Alternate  Stored  Chemical  Energy 
Reactions,"  U.S.  Department  of  Commerce  Report  NBSIR  85-3282,  1985. 

3.  PARNELL,  L.A.,  EDMUNDS,  D.G.  AND  ROGERSON,  D.J.,  "RTR  Studies  of  Closed 
Combustion  of  Liquid  Metal  Fuels,"  AIAA-87-1808,  AIAA/SAE/ASME/ASEE  23rd 
Joint  Propulsion  Conf.,  San  Diego,  CA,  1987. 

4.  PARNELL,  L.A.,  KATZ,  D.L.,  GILCHRIST,  J.T.,  BRYANT,  L.E.,  LUCERO,  J.P.  AND 
ZERWEKH,  W.D.,  "Radiography  of  Liquid  Metal  Fuel  Combustion,"  Proc.  22nd 
Intersociety  Energy  Conversion  Eng.  Conf.,  Philadelphia,  PA,  1987. 

5.  LIN,  T.F.,  MAHAFFY,  J.H.,  MANKIN,  J.C.  AND  MILLER,  T.F.,  "Thermal 
Performances  of  Oxidant  Injections  in  Submerged  Liquid  Metal  Combustion," 
27th  JANAF  Combustion  Subcommittee  Meeting,  1,  567-576,  1990. 

6.  LIN,  T.F.,  "Modelings  of  Submerged  Liquid  Metal  Combustion,"  27th  JANAF 
Combustion  Subcommittee  Meeting,  1,  577-586,  1990. 

7.  BLAKE,  T.R.  AND  PARNELL,  L.A.,  "Reacting  Gas  Jets  in  Closed  Liquid  Metal 
Combustion,"  27th  JANAF  Combustion  Subcommittee  Meeting,  1,  597-606,  1990. 

8.  CHAN,  S.H.,  JANKE,  P.J.  AND  SHEN  T.R.,  "Equilibrium  Computations  of 
Multiphase  Nonideal  Electrolytic  Systems  and  Structure  of  Turbulent 
Reacting  Dissolving  Jets,"  22nd  Symposium  (International)  on  Combustion, 
The  Combustion  Institute,  pp.  721-729,  1988. 

9.  CHAN,  S.H.,  TAN,  C.C.,  ZHAO,  Y.G.  AND  JANKE,  P.J.,  "Li-SF.  Combustion  in 
Stored  Chemical  Energy  Propulsion  System,”  23rd  Symposium  (International) 
on  Combustion,  The  Combustion  Institute,  pp.  1139-1146,  1990. 

10.  CHAN,  S.H.,  ZHAO,  Y.G.,  JANKE,  P.J.  AND  TAN,  C.C.,  "Combustion  of 
Turbulent  Gaseous  Fluorine  Jets  Submerged  in  Molten  Lithium  Fuel,” 
AIAA/ASME  Thermodynamics  and  Heat  Transfer  Conference,  Seattle,  WA,  1990. 

11.  CHAN,  S.H.  AND  SHEN,  T.R.,  "Structure  of  Turbulent  Reacting  Gas  Jets  in 
Liquids  at  Elevated  Temperatures  and  Pressures,”  Int.  J.  Heat  and  Mass 
Transfer,  Vol.  33,  No.  11,  pp.  2503-2509,  1990. 

12.  ISHII,  M.,  Thermo-Fluid  Dynamic  Theory  of  Two-Phase  Flow.  Eyrolls,  Paris, 
1975. 

13.  ABOU-ELUIL,  M.N.M.  AND  ABOU-ARAB,  T.N.,  "Prediction  of  Two-Phase  Flow  and 
Heat  Transfer  in  Vertical  Pipes,”  Proc.  5th  Symposium  on  Turbulent  Shear 
Flows,  Ithaca,  NY,  pp.  8. 1-8.7,  1985. 


-  15  - 


14.  SIRIGNANO,  W.A.,  "Heat  and  Mass  Transfer  in  Liquid-Gas  Spray  Systems," 
Presented  at  the  1990  ASME  Winter  Annual  Meeting,  Dallas,  TX,  1990. 

15.  SIMONIN,  0.  AND  VIOLLET,  P.  L.,  "Numerical  Modeling  of  Devolatization  of 
Pulverized  Coal  Injection  Inside  a  Hot  Coflowing  Air  Flow,"  Lecture  Notes 
in  Engineering  (Turbulent  Reacting  Flows),  Springer-Verlag,  NY,  Vol .  40, 
pp.  824-846,  1989. 

16.  CHAN,  S.H.  AND  ABOU-ELLAIL,  M.M.M.,  "A  Reacting  Two-Fluid  Model  for 
Turbulent  Two-Phase  Diffusion  Flames,"  ASME/JSME  Thermal  Engineering 
Proceedings,  Vol.  5,  ASME  1991,  Reno,  NV,  pp.  95-103,  March  1991. 

17.  MOSTAFA,  A. A.  AND  ELGHOBASHI,  S.E.,  "Effect  of  Liquid  Droplets  on 
Turbulence  in  a  Round  Gaseous  Jet,"  NASA  Contractor  Report  -  175063,  1986. 

18.  CHAWLA,  T.C.,  "Droplet  Size  Resulting  From  Breakup  of  Liquid  at  Gas-Liquid 
Interface  of  Liquid-Submerged  Subsonic  and  Sonic  Gas  Jets,"  Int.  J. 
Multiphase  Flow,  Vol.  2,  pp.  471-475,  1976. 

19.  ABOU-ELLAIL,  M.M.M.  AND  CHAN,  S.H.,  "Modeling  of  Turbulent  Reacting 
Multi -Phase  Flows,"  4th  Int.  Symp.  on  Transport  Phenomena  in  Heat  and  Mass 
Transfer,  Sydney,  Vol.  1,  pp.  391-401,  1991'. 

20.  CHAN,  S.H.  AND  TAN,  C.C.,  "Complex  Equilibrium  Calculation  by  Simplex  and 
Duality  Theories  with  Applications  to  Liquid  Metal  Fuel  Propulsion 
Systems,"  Combustion  and  Flame,  Vol.  88,  pp.  123-136,  1992. 

21.  PARNELL,  L.A.,  GILCHRIST,  J.T.  AND  ROGERSON,  D.J.,  "Flash  and  Real-Time 
Radiographic  Study  of  Closed  Liquid  Metal  Combustion,"  2nd  ONR  Propulsion 
Meeting,  Irvine,  CA,  pp.  188-205,  October  17-18,  1989. 

22.  WARDA,  H.A.  AND  ABOU-ELLAIL,  M.M.M. ,  "Computer  Simulation  of  Two-Phase 
Flow  in  Oil  Wells  Using  Gas-Lift,"  Engineering  Software  for 
Microcomputers,  Proc.  1st  Int.  Conf.,  Venice,  pp.  567-581,  1984. 


MASS  AND  VOID  FRACTION  •^^^SITY  (KG/M**3) 


Fi^.l  Equations  of  state  for  Li  -  SEt  system 
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Axial  variation  of  mean  centerline  mixture 
velocity  for  the  MF  and  LHF  models. 
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